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Abstract 

Neutron scattering from copper benzoate, Cu(C6D5COO)2-3D20, provides 
the first direct experimental evidence for field-dependent incommensurate low 
energy modes in a one-dimensional spin S = 1/2 antiferromagnet. Soft modes 
occur for wavevectors q = it ± Sq{H) where Sq{II) « 27rM{H)/ g^B as pre- 
dicted by Bethe ansatz and spinon descriptions of the S = 1/2 chain. Un- 
expected was a field-induced energy gap A{H) oc H", where a = 0.65(3) as 
determined from specific heat measurements. At if = 7 T {gfisH/J = 0.52), 
the magnitude of the gap varies from 0.06 — 0.3 J depending on the orientation 

of the applied field. 
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When a conventional antiferromagnet is subject to a magnetic field the antiferromag- 
netically aligned spins reorient perpendicular to the field in a so-called spin fiop transition, 
and magnetization develops through homogeneous canting of the spins in the direction of 
the field. A different scenario is expected for a spin S = 1/2 one-dimensional antiferromag- 
net where the magnetization is associated with defects in the cooperative spin-singlet 
ground state. The effect is accounted for by mapping the spin chain to a one-dimensional 
system of interacting fermions [Q, popularly referred to as spinous. The magnetic field 
causes a Zeeman splitting of the half-filled, doubly degenerate spinon band [Fig. |l](a)], and 
the spin chain develops soft modes [Fig. |l|(b)] at the incommensurate wavevectors which 
connect the field-dependent Fermi points. The new periodicity, given to lowest order by 
the ratio of the magnetic field to the spinon velocity, coincides with the average separation 
between the magnetization-carrying defects. While incommensurate spin correlations do not 
occur in classical one- dimensional spin systems such as TMMC [^, they may be a general 
feature of one-dimensional quantum spin chains [0,0. In this letter we present the first direct 
experimental evidence for incommensurate spin fiuctuations in a uniform antiferromagnetic 
S = 1/2 chain. 

To search for incommensurate spin correlations, we used inelastic neutron scattering 
to probe the wavevector-dependent spin susceptibility of the model one- dimensional S = 
1/2 antiferromagnet (AFM) copper benzoate, Cu(C6D5COO)2 ■ 3D2O. In contrast to other 
systems that have been studied in this context in copper benzoate the exchange con- 
stant J is small enough to permit the large values of the reduced field, h = qhbH/J ~ 1 
needed to observe this effect, while still being large enough to allow the interesting low 
energy part of the spectrum to be explored with a cold neutron triple-axis spectrometer 
ip. Our most important result is that a magnetic field induces new low-energy modes in 
the excitation spectrum of copper benzoate at incommensurate wavevectors g = vr ± 5g, 
where Sq{H) ^ 2'kM{H)/ gfiB, with M{H) being the magnetization per spin, as predicted 
by theory []!]-§]. The modes are not completely soft, however, because the field also induces 
a gap in the excitation spectrum both at the incommensurate wave vector and at g = tt. 
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First identified as a linear chain AFM by Date et al. [|T0[, copper benzoate is centered 
monoclinic, space group /2/c, with room temperature lattice constants a = 6.98 A, b = 34.12 
A, c = 6.30 A, and (3 = 89.5° ||rT|. Copper ions within a spin chain are separated by c/2 and 
coordinated by edge sharing, tetragonally distorted oxygen octahedra. The near-neighbor 
intrachain exchange interaction is J = 1.57 meV [^J. For this work, deuterated single crystals 
were prepared as described previously p. Specific heat measurements in magnetic fields up 
to H = 8.8 T were made on individual single crystals of typical mass 0.01 g using the 
relaxation method. For the neutron scattering measurements, approximately 500 crystals 
with total mass 3.82 g were mutually aligned to within 5° in the horizontal (hOl) scattering 
plane. Neutron scattering measurements in fields up to if = 7 T || b were performed on 
the SPINS cold neutron triple-axis spectrometer at NIST. FWHM beam divergences were 
50'/A;j(A^^) — 80' — 80' — 215', and the fixed final energy was 2.5 meV, yielding energy and 
wavevector resolution AE = 69 fieV and Ag/vr = 0.02, respectively. We refer to wave vector 
transfer along the chain as g = Q ■ c/2 = /vr. The normalized magnetic scattering intensity 
J(g, cu) 0, was derived from the detector count rate by subtracting T = 25 K data as a 
background, dividing by the squared magnetic form factor for copper, and normalizing to 
incoherent elastic scattering from vanadium. 

Figure ^ shows the g-dependence of the low energy {hu = 0.21 meV) and low temperature 
(T = 0.3 K) inelastic magnetic neutron scattering intensity from copper benzoate for four 
values oi H. At H = there is a single peak centered at g = vr. This peak is much broader 
than the instrumental resolution Ag, and arises when hut and g lie within the two-spinon 



continuum of the S = 1/2 AFM chain P,|9| jr3| , p!^ . In finite field, the data show additional 



peaks at incommensurate values of wave vector transfer, vr ± 6q, where 6q increases with 
H. The data of Fig. ^ represent the first direct experimental observation of a field-induced 
incommensurate length scale in a uniform spin-1/2 chain and is the central result of this 
paper. 

To extract the field dependence of the wavevector 6q for quantitative comparison to 
theory, one needs a specification of the full g- and ^cu-dependent magnetic density of states. 



Toward this end, we have performed both specific heat measurements, which are sensitive 
to the momentum space average of the density of states, and constant-g neutron scattering 
scans at g = TT and at g = 1.127r, which is the wavevector at which the incommensurate 
peak occurs at H = 7 T. Figure |^ shows the resulting neutron scattering data. In zero 
field at g = TT [Fig. ^(a)], the threshold for magnetic scattering is hu ~ 0, and the intensity 
decreases monotonically for hu > 0. At g = 1.127r [Fig. 0(b)], the threshold has increased 
to huj ^ 0.8 meV, with no magnetic scattering visible below that energy. The lines through 
the data in the bottom two frames of Fig. ^ were calculated by convolving the zero-field 



dynamic spin correlation function derived by Schulz [|T^ with the instrumental resolution. 

The spectrum changes dramatically for if = 7 T, as shown in Fig. ^(c) and Fig. |^(d). 
At g = vr, there is no magnetic scattering for hu < 0.1 meV, indicating that a gap has 
developed in the spectrum. Above this gap, a sharp, resolution-limited mode peaked at 
hu = 0.17 meV now marks the onset of the continuum. A second, resolution-limited mode 
also appears at g = vr at an energy close to the Zeeman energy gfisH = 0.81 meV. Figure 
^(b) reveals that at if = 7 T, the spectrum at the incommensurate wave vector g = 1.127r 
also has a gap to a resonant mode at hu = 0.22 meV. The dependence of the intensity of the 
resonant modes on the angle between the total wavevector transfer Q and the chain axis is 
consistent with the theoretical expectation [^ that the incommensurate mode is polarized 
parallel to the field, and the mode at g = vr is polarized perpendicular to it [|T^ . 

We now turn to the specific heat measurements, which explore the field dependence of 
the gap. Figure ^ shows specific heat data for copper benzoate for the same values and 
orientation of the applied field {H || la) as for the data in Fig. |^. The total specific heat for 
H = and 3.5 T is shown in the inset, plotted as Ctot/T. A small lattice contribution to 
Ctot below T = 1 K is visible in the slight temperature dependence of Ctot/T for H = 0. 
Fitting the T < 1 K zero-field data to Ctot/T = a + bT^ gives b = 0.024(5) J/mol 
This lattice term was subtracted from all the data, and the main part of Fig. ^ shows the 
resulting magnetic specific heat C. The linear magnetic specific heat we observe in zero 
field is consistent with the low-energy linear dispersion relation of the spinous in the ID 



spin-1/2 AFM. The fit shown gives C{T) = 0.68{l)R{kBT/ J), in excellent agreement with 
the theoretical value of C(T) = 0.7R{kBT/ J) ||T^. The data set an upper limit of ~ 9 fieV 
for any zero-field gap and an upper limit AS < 3 x 10~^i?ln2 on the entropy change of any 
magnetic phase transition for T > 0.1 K p!8| . 

In finite fields and at low T, C is suppressed below its zero field value. As T increases, 
C{H, T) rises above the zero- field curve, before settling back down to it at high T (inset to 
Fig. 1^). This behavior indicates a transfer of spectral weight from low to higher energies, 
and is a clear signature of a field-induced gap A{H). 

To determine the field dependence of A, we fit the H data in Fig. ^ to 

fiR f A \^'^ A , ^ n rr.^ /-,^ 



which is the low-T specific heat [|T9[ for n species of non-interacting one-dimensional bosons 



with a gap and dispersion relation 

huiq) = ^A^ + {v{q-qo)y, (2) 

A term proportional to H^/T"^ was included in the fit to take into account the small nuclear 
spin contribution. The fits are shown as solid lines in Fig. |^, and Fig. |^(b) summarizes the 
field dependence of the gaps derived from this analysis as well as for data we have taken 



with H along the other principal magnetic directions, c" and a" |]T2[. The average of the 
gaps measured by neutron scattering for = 7 T || b, at g = vr and 1.127r (filled symbols in 
Fig. 0(b)) correspond nicely to the gap derived from specific heat. Moreover the values of 
v/n = 0.49(3), 0.55(3), and 0.96(5) meV/mode for H \\ b, c" and a", respectively, are not 
far from the crude estimate v/n ^ 0.57rJ/6 = 0.41 meV/mode, where n = 6 is the number 
of soft modes in a Brillouin zone. The field dependence of the gaps is described by the 
power-law A{H) = AH" with a = 0.65(3). The prefactors A are in ratios 1 : 2.0 : 0.55 for 
H applied along the b, c" and a" directions, respectively. A slightly better fit to the c" data 
is obtained by including a finite critical field: A{H) = A{H — i/c)", where He = 0.25(3) T 
and a = 0.58(3). We note that small field-dependent gaps are not found in the classical easy 



plane AFM spin chain TMMC, where a field in the easy plane induces a central diffusive 
soliton mode and a gap A{H) = gjisH PJ20|]. 

With the information about the field-dependence of the gap derived from the specific 
heat data, we can extract the field-dependent incommensurate wave vector from the neutron 
scattering data in Fig. 0. We approximate the low energy response as resonant modes with 
dispersion relations given by Eq. (||), where go = or tt ± Sq{H) (see Fig. |l](b)). The 
velocity v was fixed at the zero field value 7rJ/2 |^T|, and the gap A{H) was taken from 
the specific heat measurements. Although it neglects the continuum scattering above the 
resonant modes that is visible in Figs, ^(a) and 0(b), this model provides a reasonable fit 
for the finite field const ant-Zicj scans in Fig. |^. Note that for if = 3.5 T, it appears that 
a single branch of the dispersion relation at the incommensurate wave vector is resolved, 
while the other merges into a broad central peak at g = vr. Thus the incommensurate 
wavevector at this field actually lies at the local minimum between the peaks. This is 
in contrast to the H = 5 T and H = 7 T data, where the incommensurate modes are 
well-separated from the q = it mode, and the energy of the scan is very close to the gap. 
The fits establish that the incommensurate peaks are resolution-limited, and also provide 
a systematic way of deriving the field-dependent incommensurate wavevector 6q{H). This 
is shown as a function of H in Fig. |^(a). The solid line in this figure is the theoretical 
prediction Q: 5q{H) = 2t:M{H) / g^iiB with g^, = 2.06 0. M{H) is calculated using Eq. 



(2.24) of Ref. 0, which provides an excellent description of the magnetization for copper 
benzoate ||22|. Despite the complications associated with the field-induced gaps, the theory 
accounts well for the field-dependent positions of the incommensurate soft modes. 

In summary we have provided the first direct experimental evidence for field-dependent 
incommensurate soft modes in a S=l/2 one dimensional antiferromagnet. Thus, one of the 
key features of the spinon description of S=l/2 AFM chains, namely a chemical poten- 
tial which can be regulated by an external field, has finally been verified. Our results also 
have implications beyond the nearest-neighbor coupled S=l/2 chain. For example, theoret- 



ical studies predict that both long-range coupled models p3| and systems such as ladders 



and dimerized chains with intrinsic gaps should show similar behavior, indicating that 
field-induced incommensurabilities are a general property of quantum spin chains. The in- 



commensurate lattice distortions recently observed in TTF-CuBDT and CuGeOs p5 
may therefore be closely related to our results. In copper benzoate, the field also introduces 
gaps in the excitation spectrum both at the commensurate and at the incommensurate wave 
vector, something which was not predicted by theories for isotropic Heisenberg spin chains. 
Possible causes for the gaps are diagonal as well as off-diagonal exchange anisotropics, and 
a staggered g-factor anisotropy [|l2l which induces a small effective staggered field when a 
homogeneous external field is applied . 
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FIGURES 

FIG. 1. (a) Schematic of the Zeeman-spht bands of the fermion mean field theory for the 
S=l/2 AFM chain in a magnetic field H. The Fermi energy is at = 0. Solid (dashed) arrows 
show incommensurate (commensurate) spanning vectors, (b) Boundaries of the resulting spinon 
continuua, for spin fluctuations parallel (solid lines) and perpendicular (dashed lines) to H. Note 
incommensurate soft modes near q = n. Horizontal dotted line shows trajectory of scans in Fig. 2. 

FIG. 2. Magnetic scattering at huj = 0.21 meV along the (0.3,0,/) directionfor four values of 
magnetic field H at T = 0.3 K. The solid lines for H ^ are fits to resolution-limited spin wave 
scattering as described in the text. The line through the H = T data is the theory ||l^ for the 
dynamic spin correlation function convolved with the experimental resolution. A correction factor 
varying from 0.7 to 0.95 over the range of the scan was applied to the model to take into account 
the measured q dependent sample transmission. 

FIG. 3. Energy dependence of the magnetic scattering intensity at T = 0.3 K for Q = (0.3, 0, 1) 
and Q = (0.3,0,1.12) at T = 0.3 K. The latter wave vector corresponds to the position of the 
incommensurate maxima in the H = 7 T constant-/iw = 0.21 meV scan. The solid lines through the 
H = OT data are the theoretical dynamic correlation function ||l^ convolved with the experimental 
resolution. The dashed lines are the average of the data in (b) for 0.08 meV < huj < 0.6 meV, 
which is a good measure of the over-subtraction caused by isotropic magnetic scattering contained 
in the T = 25 K data used as a background. 

FIG. 4. Specific heat of copper benzoate as a function of temperature for four values of applied 
field H II h. Inset: total specific heat plotted as Ctot/T for H = and 3.5 T. Main panel: magnetic 
heat capacity C after subtraction of phonon contribution. Solid lines are from fits described in the 
text. 
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FIG. 5. (a) Field dependence of the displacement 5q of the incommensurate side peaks from 
g = TT in copper benzoate, as derived from fits to the data shown in Fig. 2. The solid line is 
the theoretical curve from Ref. Q. (b) Field dependence of the energy gap derived from fits to 
specific heat data such as those shown in Fig. 4. Data for fields along the three principal magnetic 
directions are shown. Filled symbols are the gaps measured at q = n and q = 1.127r by neutron 
scattering. The solid lines are from fits to power-laws described in the text. 
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Fig. 1 Dender et al. 
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Figure 2. Dender et al. 
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Fig. 4, Dender, et al. 
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Fig. 5 Dender, et al. 



